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ABSTRACT 

A recently discovered filament of polarized starlight that traces a coherent 
magnetic field is shown to have several properties that are consistent with an ori¬ 
gin in the outer heliosheath of the heliosphere: (1) The magnetic field that pro¬ 
vides the best fit to the polarization position angles is directed toward £=357.3°, 
b= 17.0° (±11.2°); this direction is within 6.7° ± 11.2° of the observed upwind 
direction of the flow of interstellar neutral helium gas through the heliosphere. 

(2) The magnetic field is ordered; the component of the variation of the polar¬ 
ization position angles that can be attributed to magnetic turbulence is ±9.6°. 

(3) The axis of the elongated filament can be approximated by a line that defines 
an angle of 80° ± 14° with the plane that is formed by the interstellar magnetic 
held vector and the vector of the inflowing neutral gas (the BV plane). We 
propose that this polarization feature arises from aligned interstellar dust grains 
in the outer heliosheath where the interstellar plasma and magnetic held are de¬ 
flected around the heliosphere. This interpretation suggests that the polarization 
is seen where stream lines of the how and the draped interstellar magnetic held 
lines are approximately parallel to each other and perpendicular to the sightline. 

An ordered magnetic held is required so that grain alignment is not disrupted 
during the interaction. The filament location is consistent with dust plumes 
previously predicted to form around the heliosphere. The proposed outer he¬ 
liosheath location of the polarizing grains can be tested with three-dimensional 
models that track torques on asymmetric dust grains as they propagate through 
the hcliosheath plasma, and using these models to evaluate grain alignment and 
the asymmetric extinction of the grains. 

Key words: Sun: heliosphere - ISM: magnetic helds - ISM:dust - polarization 
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INTRODUCTION 


The Sun’s heliosphere and the astrospheres of other stars sweep up interstellar gas and 
dust as they travel through interstellar clouds. Infrar ed data s how that dusty bow waves, 
bow shocks, and shells are common around other stars flPeri et al.l 2012b Three-dimensional 
MHD models of the interactions between interstellar dust grains and the heliosphere predict 
that dust plumes com posed o f sub-m icron grains with high charge-to-mass ratios form aroimd 
the heliosphere flSlavin et al.l 2012). Detection of such a dusty shell around the heliosphere 
using infrared measurements would be challenging because of the difficulty in separating 
background emissions from a heliospheric feature. 

The center of the “ribbon” arc of energetic neutral atoms (ENAs), discovered by the 


ing over the heliosphere ( 

McComas et al. 

2009; 

Schwadron et al. 

2009; 

Funsten et al. 

2009, 

2013. 

2014; 

Heerikhuisen et al. 

2014 

). The extraordinary circularity and symmetry of the 


ribbon supports a hypothesis that the interaction between the heliosphere and the inter¬ 
stellar medium (ISM) involves an ordered magnetic field. Charged interstellar dust grains 
interacting with the ordered field may polarize starlight and provide an alternate method 
for sampling a dusty shell around the heliosphere. 

The Sun is traveling through an interstellar medium (ISM) con taining part ially ionized 
interstellar gas with low average densities, a magnetic field, and dust (IFrisch et al.l 2 0111) . The 
electron density, ~ 0.03 — 0.1 cm -3 , and magnetic field strength, ~ 3 /iG, in the surrounding 
interstellar gas are typical of the mag netoionic medium that contains a stru ctured magnetic 

About 0.5%- 


field and fills 10% — 40% of space flHaverkornl 12010 ; H affner et al 


2009) 


0 .7% o f the local interstellar mass density is in charged dust grains (ISlavin fe Frisch! 2008 
K ru ger et al. 2014). Mapping the structure of the nearby magnetic field requires high 


sensitivity measurements of starlight that is polarized by asymmetric dust grains aligned 
with respect to the foreground interstellar magnetic field (ISMF). 

This is the fourth study in a global survey of the polarizations of light from nearby stars 
in order to map the ISMF w ithin 40 pc an d compare it with the direction of the ISMF shaping 
the heliosphere (IFrisch et al.l 2010, 2012 . 12015 . Papers I, II, III). The survey includes stars 
within ~ 90° of the galactic center. A magnetic filament was identified in these polarization 
data through the linear rotation of their polarization position angles with distance (Paper 
III). The filament forms a band over 100° long and ~ 15° wide in the sky, and extends to 
within 10 pc of the Sun. Paper III discusses a possible association of this filament with 
the Loop I superbubble. Evidence that this magnetic filament may be associated with the 
heliosphere is the topic of this paper. 
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In the diffuse interstellar medium where the ISMF direction is also traced by the Faraday 
rotation of radio waves and the polarization of synchrotron emission, the optically asym¬ 
metric interstellar dust grains that linearly polarize starlight are aligned so that plane of 
polarization is parallel to the ISMF direction, and is strongest where B • R=0 for magnetic 


field B and radial sightline R (e.g. iDavis &; Bergd Il968l; iLazarianl 120071; lAndersson et al. 
2014 1. The polarization position angled of the survey data serve as a method for evaluating 


ISMF direction that is traced by the polarization the data. In Paper III, the filament star 
polarizations were found to track a separate ISMF direction than the other polarization data 
in the sample. After the hlamanet subset of data was excluded from the analysis of the ISMF 
direction, the dominant galactic magnetic held in the immediate solar vicini ty was f ound to 
have a magnetic pole located at £ = 36°, b = 49°, (±16°) (Paper III, Frisch et ah 2014 1. 
This dominant local ISMF direction agrees to within 1-sigma with the ISMF d irecti on trace d 
by the weighted-mean center of the highly circular IBEX ribbon arc (jFunsten et ah 2013)). 
The nominal pole directions differ by 7.6° ± 16.2°. 

The velocity vectors of interstellar dust and gas fl owing into the inner heli osphere are 
the same. Dust a rrives from I = 8°, b = 14° f ±13°l (IFrisch et al.l Il999l: iBelheouane et al. 
2012 ; Mann 2010 1 at a velocity of -24.5 km s _1 

6=16°, at a velocity of 


2004; 

Mobius e 

al. 

■_ ~ ?. 

2012 

McComas et al. 

2014b|). 


26 km s 1 (Bzowski et al. 2014 iWitte 


Bzowski et al.l 120121; IWood et al.l I20l4 iKatushkina et a, hi 12014 ; 


flKimura et al.1 120034 . Neutral helium arrives 


There is a deficit of low-mass interstellar grains (< 10 -13 gr) in 
the inner heliosphere compared to the parent interstellar population that is attributed to 
the action of Lorentz forc es on grains with high charge-to-mass ratios in the outer heliosphere 


and heliosheath regions flKimura X. Mannl Il998l; IFrisch et al.l Il999l ; iLinde fe Gombosil 12000 


Slavin et al.ll2012l;lSterken et al.ll2012l ) . Three-dimensional models of the trajectories of inter¬ 


stellar dust grains interacting with a three-dimensional MHD heliosphere model predict that 


t hese excluded grains create shells or plumes of dust around the heliosphere (iSlavin et al. 

2012h . 


2. PROPERTIES OF THE MAGNETIC FILAMENT 


An ordered component of the local ISMF within 40 pc was originally noticed in Paper II, 
using the high-sensitivity polarization data set acquired with PlanetPol (iBailev et al.l 20 101). 
The distinguishing characteristic of this ordered ISMF was that the polarization position 


lr The polarization position angle is the angle between the plane of polarization of the E-vector and the 
north pole, and labeled positive for rotations of the postion-angle meridian toward larger longitudes. 
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angle in the equatorial coordinate system (# C ei) rotated linearly with the distance of the 
star. Using an expanded set of stars that included the PlanetPol data, we found a linear 
relation between the polarization position angle in galactic coordinates, 9 ga i, and the star 
distance (Paper 111). The stars separated into two groups that had slightly different distance 
dependences for 9 ga i, neither of which agreed with the gradient found from 9 ce \ in Paper If. 
This difference was a clue that the apparent relation between polarization position angle and 
distance was spurious, and that another factor organized these polarization position angles. 
The stars showing the regular variation of position angle, $pa, with distance were found to 
trace an extended filamentary feature of linearly polarized starlight that extends over 100° 
in the sky. The polarization data tracing this filament are listed in Table [1] and plotted in 
Fig. [T| in a nose-center ecliptic coordinate system^. 

In order to search for the organizing principle of the filament, data on thre e new stars in 


the region of the filament were collected with DiPol-2 at the UH88 telescope flPiirola et al. 


2 014 1. These new data are also listed in Table 1. These polarization data that define the 
filament were acquired using different instruments at seven observatories; the mean errors of 
the measurements are not uniform between these different data sets. This extended set of 
polarization data of sixteen stars is used here to search for the controlling factor that creates 
the alignment of the filament polarizations. Several stars have multiple measurements; all 
data were used in this analysis. This selection of data is relatively unbiased in the sense that 
only the three stars observed most recently by DiPol-2, HD 153631, HD 141937, and HD 
145518, were observed for the purpose of confirming the filament properties. 


The first test performed on the filament polarizations is to determine the direction of the 
ISMF associated with these polarizations. The best-fitting 1SMF to the polarization data 
in Table [T] was obtained with the method in Papers If and 111, where it is assumed that the 
ISMF direction can be described by a great circle in the sky, and that the best-fitting ISMF 
direction is the one that minimizes the ensemble of sines of the polarization position angles 
referenced to the ISMF great circle^ A polarization position angle that is perfectly aligned 
with the ISMF direction will have a position angle of either 0° or 180° when referenced to 
the ISMF direction, so that the sine of the polarization position angle, 0 ism f, is zero. The 
polarization strengths are not used as the basis for determining the ISMF direction. 


The heliosphere no se is defined by the dire ction of inflowing neutral interstellar helium atoms, A=255° 


/ 3 = 5°, or £= 4 °, 6=15° iMcComas et al.ll2014bh 


3 In the context of the software, this assumption is implemented by rotating the coordinate system into 
a 1-degree grid covering the sky to find the direction that best-matches the rotated polarization position 
angles. No assumption is made about the three-d imension al di stribu tion of the ISMF since the polarizing 
dust screens can be anywhere in the sightline. See [Frisch et all ( 201 5) for more details. 
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Testing this set of position angles against all possible ISMF directions, using a 1-degree 
grid and the weighted means of the position angles to evaluate the minimum, yields a best¬ 
fitting ISMF direction to the polarizations of the filament stars that is toward £=357.3°, 
6=17.0° (±11.2°). The uncertainties on this direction are obtained from the standard devi¬ 
ation of the polarization position angles calculated with respect to the magnetic held frame, 
6*i sm f, that corresponds to the best-fitting magnetic held pole. An estimate of the magnetic 
turbulence of this ordered held is obtained by assuming that the uncertainty of 11.2° is com¬ 
posed of the root-mean-square of the turbulent held component and the mean uncertainty 
of the position angle measurements, A$pa=5.6°, so that the “turbulent” component of the 
position angles is approximately ±9.6°. 


The surprising result of best-fitting ISMF direction to the filament polarizations is that 
it is within 6.7°±11.2° of the upwind velocity vector of the primary interstellar neutral helium 
how into the heliosphere, whic h arri ves from the direction of £=3°, 6=16° (or A=255°, (3= 5° 
in ecliptic coordinates, e.g. IMcComas et al. 2014q ). The coincidence between the directions 
of ISMF vector and the velocity vector of neutral interstellar gas flowing into the heliosphere 
suggests that the filament may be related to the heliosphere, and that the hlament axis may 
trace a feature of the heliosphere. 


The major axis of the hlament axis can be approximated by fitting a line through the 
galactic coordinates of the stars. The best-fitting line corresponds to £ = 25.118° (±1.258°) + 
0.139° (±0.037°) x S£, with a reduced y 2 of 1.14; here S£ is the difference between the star 
galactic longitude and the end point of the polarization hlament at £=67° (Table 1). The 
long axis of the hlament is shown as the gray line in Fig. [2j 


The angle between the hlament axis and the BV plane of the heliosphere is 80° ± 14°, 
so that the hlament axis and BV plane are perpendicular. The BV plane, projected as the 
gray dashed line in Fig. [21 is the plane formed by the ISMF (B) and neu t ral interstellar win d 
velocity (V) vectors that shape the heliosphere flPogorelov et al 2QQ9; Opher et al, 2009'). 
For the magne tic held direction (B) we use the center of the IBEX ribbon arc at 4.3 keV 
fiFunsten et al.l 120131 ). At this energy ENAs will have longer mean free paths than those 
at lower energies, so that this direction extends further into the ISMF than the directions 
from lower energies. The interstellar velocity (V) i s given by t he velocity vector of neutral 
interstellar helium flowing through the heliosphere flMcComas et al.li2014bj ). Ribbon models 
indicate t hat the ri bbon cent er is within 5° of the true ISMF direction for likely held strengths 
of 3 /iG f Heerikhuisen et al. 2014). 


A range of heliospheric data indicate that the BV plane plays a fundamental role in 
organizing the interaction between the ISM and the heliosphere. Charge exchange coupling 
between interstellar neutrals and ions cause an offset of the neutral gas flow along a direction 
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that is parallel to the BV plane; examples are the 4.7° deflect ion the interstellar hydrogen 
flow from the benchmark He° flow direction flLallement et a 
secondary helium population (a "warm breeze”, Ku biak et ah 


20101) . and the offset of the 


2014 . Fig. [2|). The BV plane 


also pl ays an organi zing role in the spectral dependence of the ribbon symmetry, with ENAs 
(iFunsten et al.l 2014 h and organizes the heliosphere asymmetries (Fig. [3]). 


Figure [Q plots the filament polarizations over a simul ation of t he co lumn densities of 
0.01-1.0 /ini interstellar dust grains within 400 AU from lSlavin et al.l 02012 . S12). The figure 
is an ecliptic projection centered o n the_nose of the heliosphere, defined by the inflowing 
interstellar helium gas. ISlavin et al.l (120121 S12) modeled the tr ajectories of interstella r 


compact silicate grains through the 3D M HD heliosphere models of iPogorelov et al.l (120091) . 
These models were successfully used by Sc hw adr on et al. (j,2009[) to show that the IBEX 
ribbon appears in sightlines that are perpendicular to the ISMF draping over the heliosphere, 
i.e. where B • R=0 for magnetic field vector B and radial direction R. The models include an 
ISMF direction that is within ~ 17° of the interstellar field direction traced by the weighted 
center of the IBEX ribbon arc@ Densities of 0.01 /.irn interstellar grains within ~ 400 AU of 
the Sun are shown in Figure Q] based on these models. The filament crosses the heliosphere 
nose at the plot center and reaches high latitudes on the port side of the heliosphere at HD 
172167, and low latitudes on the starboard side of the heliosphere at HD 119756 (borrowing 
nautical terminolog}@ ) The BV plane controls the symmetry of the cloud of interstellar 
dust grains interacting with the heliosphere, and the asymmetry of the polarization filament 
about the BV plane is evident. The filament is slightly more twisted toward the vertical 
direction in the north and port side of the heliosphere compared to the starboard side in this 
simulation. 


4 The gala ctic coordinates of the weighted mean center of the IBEX ribbon arc are i = 36.7° ± 2.1°, b = 
56.0° ± 0.6° ( Funsten et al. 2013 : Frisch ct al. 2014 ). 


5 In nautical usage, port and starboard are referenced to the bow of a ship. Using an xyz right-handed 
cartesian coordinate system, the x-axis is directed towards the bow, the y-axis is directed toward the port 
side and the z-axis points ”up” toward the opposite direction of the gravitational force. By analogy, for the 
heliosphere the nose-direction of the heliosphere corresponds to the x-axis, the y-axis is directed toward the 
port side of the heliosphere, and the z-axis, pointing north, completes the right-handed coordinate system. 



































3. DISCUSSION 


3.1. Possible Origin of Filament Polarizations in Outer Heliosheath 


The agreement between the ISMF direction of the filament polarizations and the velocity 
vector of interstellar material flowing through the heliosphere, the perpendicularity between 
the BV plane and the filament axis, and the low level of magnetic turbulence in the filament, 
suggest that the filament polarizations may be associated with the deflection of charged 
grains, the ISMF, and interstellar plasma around the heliosphere in the outer heliosheath 
regions. 


Interstellar nanograins with large charge-to-mass ratios ( Q/m ) ar e deflected i n the 


outer and inner heliosheath regions, and throughout the heliosphere ( 

Frisch et al. 

1999; 

Kimura & Mann 

1999 

: Linde & Gombosi 

2000; 

Kruger et al. 

2014; 

Czechowski & Mann 

2003 

Mann & Czechowski 

2004; 

Slavin et al. 

2012|; 

Sterken et al. 

2011 

2). Most grain trajectory 


simulations have assumed compact grains. Tiny grains are subject to the small parti cle 
effect , where the efficiency of secondary ejection of electrons is enhanced (iKimura fe Mann 


19981) . Porous or fluffy grains that consist of aggregates of tiny grains are also subject to 


this effect, which leads to higher grain charges and increases the a verage mass of the outer 
heliosheath grains above the predicted levels for compact grains flMa et al.l 20131). These 
deflected grains are candidates for the polarizing dust grains in the filament. 


The mass distributions of the interstellar dust grains observed w i thin 5 AU of the Sun 
exhibit a deficit of low mass grains, < 10~ 13 gr (IFrisch et al.l 11999c lLandgraf et al.l 12000 


Kruger et al 


Mathis et. al. 


20141) . when compared to the nominal power-law interstellar distribution of 
fll977l ). The deficit of small grains is attributed partly to the “filtration” of 
grains with large values oiQ/m th at cou ple to the int erste llar plasma and m agne tic field that 


is sw e pt around the heliopause (IFrisch et al 


2000; iKriiger et al.ll2014 ; 

20121 ). 


1999 


Kimura fc Mann 


1999 


Czechowski fc Mann 2003^ Mann fc Czechowski 


Linde fe Gombosi 


2004; Slavin et al. 


The perpendicular relation between the filament axis and the BV plane would be an 
outcome of the Lorentz force on the grains. Nanograins experience additional charging in 
the outer heliosheath due to collisions with the heated compressed interstellar plasma and 
enhanced radiation field, with the maximum differential increase in the surface potential of 
the grains occurring for the smallest grains (Fig. 2 of S12). S12 showed that the surface 
potentials of compact interstellar silicate grains in heliosheath regions are enhanced by an 
order of magnitude compared to values in the surrounding ISM. As the dust grains approach 
the heliosphere from the distant ISM, Lorentz forces act to displace grains into directions 
that are perpendicular to the BV plane. The action of the Lorentz force on grains with large 
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Q/m the grains naturally explains the orientation of the filament axis with respect to the 
BV plane. Fig. [3] shows the gyroradius of a 0.01 //in radius interstellar grain propagating 
through the outer heliosheath region, and as the grain is deflected around the heliopause by 
the Lorentz force; the grain trajectory model is based on the simulations in S12. 


The ~ 15° width of the filament is not yet explained. Since all of the smallest grains are 
deflected, the filament width must depend on the detailed configuration of the heliopause, 
the grain charging that occurs in the outer heliosheath, and the strength of the Lorentz 
force. The dust plumes in the models of S12 (Fig. 13 of S12, left) have a limited width that 
is similar to the filament width, and a similar spatial orientation perpendicular to the BV 
plane. Those dust plumes therefore provide a plausible model for the polarization filament 
location, provided that the efficiency of polarization can be explained. 


If the filament has a heliospheric origin, then the second characteristic of the filament, 
that the ISMF is directed toward the inflow vector of interstellar neutrals, implies that the 
filament forms where the interstellar ion and dust streamlines are parallel to each other and 
to the localized ISMF direction. The surroundi ng interstellar gas consists of about 23% ions 
and 77% neutral atoms (ISlavin fe Frischl 120081) . The collisional mean-free-path in the ISM 


approaching the heliosheath regions is r^j 330 AU flSpangler et al l 


thickness of the plasma density ramp outside of the heliopause ( Zank et al. 20131) . @ 


2011), and is similar to the 


Interstellar ions and neutrals follow different trajectories in the outer heliosheath. Ions 
are deflected around the heliopause, and the neutrals propagate through the plasma, with 
some loss due to the formation of ENAs. The frozen-in interstellar held becomes fully 
deflected where the ions decouple from the inflowing interstellar neutrals at the heliopause, 
so that the deflected plasma streamlines and magnetic held should be directed toward the 
undehected neutral velocity vector. A correspondence between the ISMF direction and the 
He° how velocity suggests that the directi on of t he drape d ISMF has twisted into the direction 
of the dehected gas how (e.g. see lOp he r fe Drake 2013). 


An important requirement is that the dust grains remain aligned with respect to the 
ISMF direction during the deflection of the dust how at the heliopause. The organized set 
of filament polarizations must be tracing an organized magnetic held that has a minimal 
turbulence. The recent in situ measurements of the magnetic held in the inner and outer 
heliosheath region s with instruments on board Voyager 1 support an ordered h eld in the 
outer heliosheath (IBurlaga et al.ll2013t iBurlaga fe Nessll20l4 iBurlaga et al.lI2014T) . Voyager 
1 crossed into the interstellar magnetic held in mid-2012 and found a more ordered interstellar 


6 The inn er heliosheath thickness found from Voyager 1 and IBEX data is 
Hsich ct al. 20ldh . 


27 AU ( Gurnett ct al. 12013 : 
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magnetic field in the outer heliosheath than had been seen in the inner heliosheath, with less 
turbulence, and a maximum standard deviation of the magnetic held strength that was 3% of 
the total held strength. In contrast, in the inner heliosheath the maximum standard deviation 
of the compressive turbulence (corresponding to variations in magnetic held strength) was 
36% of the total held strength suggesting that the inner heliosheath is unlikely to be the 
origin of the polarizations. 


The low level of interstellar held turbulence in the outer heliosheath found by Voyager 1 
suggests that high magnetic pressures will concentrate the nanograins in the outer heliosheath 
and increase the efficiency of the grains to linearly polarize light. However, Rayleigh-Taylor 
instabilities can disrupt the helio pause and allow the charged interstellar part i cles t o pen¬ 
etrate into the inner heliosheath flBorovikov fc Pogorelovl 12014 lAvinash et al .11201 41). The 


blunt heliosphere inferr ed from Voyager data of th e neutral current sheet fjBurlaga fe Ness 


2014) and plasma flow flRichardson fe Decker! 120141). and t h e exten ded high-pressure region 


at the heliosph ere nose seen in IBEX data flMcComas et al.112014al) and predicted by helio¬ 
sphere models (IPogorelov et al.l 2 0111 ). is going to affect the draping of the ISMF over the 
heliosphere in ways that are not yet fully understood. 


A ’cartoon’ illustrating the geometric relation between the filament polarizations and 
the heliosheath is^ shown in Fig. [3] the surface of the heliopause is from the MHD model of 


Pogorelov et al . (2 01 lL Fig. 1), where the BV plane and Voyager trajectories are also shown. 


The filament orientation perpendicular to the BV plane is displayed, with the filament 
located close to the heliopause on the upwind side. The hypothesized local distortion of the 
magnetic held direction parallel to the filament and heliopause is also illustrated. 


3.2. Alignment of the Filament with respect to the IBEX Ribbon 

For starlight that is linearly polarized in the ISM, a polarization position angle that 
is aligned with the ISMF direction will yield the strongest polarizations where the ISMF 
is perpendicular to the radial sightline, e.g. where B • R=0. In contrast, the filament 
polarizations trace a magnetic held direction that is directed toward the upwind velocity 
of neutral interstellar gas instead of the ISMF direction corresponding to the center of the 
IBEX ribbon arc, ~ 48° away from the gas velocity. Our hypothesis is that the polarization 
filament is visible where the dust grains are entrained in the ISMF that has been dehected 
around the heliosphere, so that polarizations would be strongest where the held lines are 
quasi-perpendicular to the sightline for otherwise-equal conditions. 

Comparison between MHD heliosphere models and the IBEX ribbon location show that 
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the ~ 1 keV ribbon is observed in sightlines where the interstellar magnetic field becomes 
perpendicular to the sightline as the field drapes over the heliosphere, i.e. where B • R=0 


flSchwadron et al. 


2009 


Schwadron fc McComas 


: VlcComas et al.l 120091: iHeerikhuisen et al.l l2010l; IChalov et al.l 12010 


2013]). Fig. [2] shows that the filament polarizations are roughly 


parallel to the IBEX ribbon. 

There is no consensus on the formation scenario for the IBEX ribbon, although the orig¬ 
inal model where the ribbon is viewed along sightlines perpen dicul ar to the_ISMF draping 
over the heliosphere is generally accepted (for a review see lMcComas et al.l 2014c). Formation 
from secondary ENA’s upstream of the heliopause is currently the leading contender for ex¬ 
plaining the ribbon. In these models the ribbon ENAs form upstream of the heliopause where 
charge-exchange between interstellar neutral hydrogen atoms and the pickup-ion plasma, 
which came from the neutralized outward moving solar wind, create energetic neutral atoms 
(ENAs). These ENAs propagate back toward IBEX if the ion pitch angles retain the domi¬ 
nant radial momentum of the original outward flowing particles. Models relying on secondary 
pi ckup i ons to generate t he ribbon place the origin 50 — 100 AU upstream of the heliopause 
(iHeerikhuisen fe Pogorelov 2011 1 up to several hundred AU. Presumably the interstellar ions 
and the dust grains with small Q/m will couple more tightly to the deflected ISMF than 
do the pickup ions that create the neutral ENAs, so that the polarization filament would 
appear at the region where the small Q/m grains are fully deflected rather than precisely 
where B • R=0. 


3.3. Are Outer Heliosphere Grains enough to Polarize Starlight? 


Interstellar column densities in the outer heliosheath are expected to be lower than total 
column densities toward the nearest stars. The gas column densities in the outer heliosheath 
are only 0.06%-0.9% of the total interstellar column densities toward t he nearb y upwind 
star 36 Oph (6 pc), if the hydrogen-wall and interstellar data in Wood et al. (120001 ) are used 
to benchmark this relation. For an average h ydrogen- wall density of n(H°)=0.25 cm -3 , the 
three-component hydrogen-wall model of Wood et al. (120001) suggests an outer heliosheath 
thickness of ~ 100 AU toward 36 Oph. If the filament polarizations are formed in an outer 
heliosheath region that is on the order of 100 AU thick, then either grain alignment is 
significantly more efficient in the heliosheath than in the local ISM, or an unknown factor 
such as foreground depolarization is reducing the interstellar polarizations of nearby stars. 


The outer heliosheath physical conditions differ from those in the ISM upstream of 
the heliosphere. Although detailed models are required to evaluate the alignment of grains 
propagating through the outer heliosheath, the factors that increase polarization efficiency 
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in this relatively thin region can be described qualitatively. These factors include higher 
grain surface potentials, enhanced strength and asymmetry of the radiation field, enhanced 
magnetic field strengths, low levels of magnetic turbulence, and the absence of foreground 
depolarizations, all when compared to the values in the local ISM. The presence of fluffy or 
porous grains, which are more sensitive to secondary electron ejection than compact grains, 
further enhancing grain alignment. 


Self-consistent simulations of compact interstellar grains interacting with the heliosheath 
regions show that surface potentials of the grains increase by up to an order of magnitude in 
the outer heliosheath, depending on grain radius (S12). The higher Q/m values will increase 
grain charge and the Lorentz force on grains, and minimize possible collisional disruption of 
the grain alignment. 


Rad iative torques p la y a ro le in grain alignment in interstellar clouds flAndersson et ah 


2014). Hoang fe Lazarianl (12014T ) have modeled the polarizations expected for the very local 
ISM, taking into account that radiative torques on optically asymmetric grains play a domi¬ 
nant role in grain alignment in the warm local ISM where gas column densities are typically 


10 1 ' ' — 10 18 '' cm -2 . Their models were able to account for polarization strengths as low 
as 0.002%, depending on the total column density of the gas. The polarization strengths of 
the filament stars range from 0.0032% to 0.1% (Table [1]), however the filament polarizations 
probably form over a column that is < 0.002 times the length of the column of interstellar 
polarizations. Radiative torques that enhance grain alignment will become of increasing 
importance as grains approach the heliopause and experience a stronger, more asymmetric, 
and more energetic radiation field (see Fig. 3 in S12). 

As the ISMF strength increases as field is compressed against the heliopause by the solar 
motion through the surrounding ISM, interstellar dust grains will couple more tightly to the 
ISMF because of smaller gyroradii. The interstellar field strength is enha nced b y factors of 
three or more outside of the heliopause in the models used byl Sla vin et al. ( 20121 ). Voyager 1 
found that field strengths in the outer heliosphere reach 6 /iG, which suggests a compression 
factor of tw o when compared to the 3 uG field strength of the interstellar cloud around the 
heliosphere fjBurlaga fe Nessl I20l4 ISlavin fe Frischl 120081: ISchwadron et al.l 120111) . 


Polar izations are not additive in the ISM since f oreground clouds may depolarize polar¬ 
ized light ( Serkowski et al.lll975 ; Davis &; Bergelll968 ). There is no information on foreground 
depolarization in the local ISM so the local significance of this effect is unknown. However, 
when compared to interstellar values, there is unlikely to be foregroimd depolarization of 
polarizations originating in the outer heliosphere. In principle the absence of depolarization 
will further increase the efficiency of heliosheath polarizations above those in the ISM. 
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A significant enhancement of polarizations compared to previous simulations of polar¬ 
ization efficiencies may arise from the physical structure of the excluded grains, and whether 
they are compact or porous aggregates (“fluffy”). Photoionization models of interstellar gas 
around the heliosphere show the presence of olivine silicates and the absence of carbonaceous 
grains (jSlavin fe Frisch] 120081 ). The grain trajectory models that predict dust plumes around 
the heliosphere (S12) are based on compact olivine silicate grains with densities of 3.3 gr 
cm -3 . In situ measureme nts by Stardust found i nterste ll ar gra ins with very lo w vo lume 
densit ies, < 0.07 gr cm 


-3 


WestDhal & the Stardust Team 

2014 

). Recently, 

Sterken et al. 


(120151 1 found an indication of porous grains from comparisons of the flow directions of in¬ 
terstellar grains in 16 years of interstella r dus t measurements from Ulysses , such as was 
proposed in earlier discussions (cf. Marini 120101 ). In situ measurements of interstellar dust 
in the inner hel iosphere have found large grains with masses 2 x 1Q~ 10 gr (radii ~ 2 /jm for 
compact grains, iFrisch et al.lll999l: ISlavin et al.ll2012t iKriiger et al.ll2014r) . These heliospheric 
large grains are similar in size to the lar ge micron-sized aggregate grains that are observed 
in the cores of dense interstellar clouds (IPagani et al. 2010; ISantos et al. 2014). Outside of 
dense clouds the l arge grain s couple to the interstellar gas over pathlcngths of ~ 50 — 100 
pc (IGrun fe Landgrafl 200 01). and the local population of large grains may be enriched by 
porous dust grains from fragmented dense clouds. 

Porous grains have larger cross-sections than the canonical compact grains that have 
been used for modeling the dust-heliosphere interactions. iMa et al. (120131 ) have modeled the 
grain charging of aggregate grains and found that the surface potentials of porous grains 
is increased by 30% or more in the outer heliosheath, which allows a stronger coupling 
between the dust and deflected heliosheath plasma and magnetic held. Aggregate porous 
grains have higher charges than compact grains, and so are deflected differently. The high 
end of the cutoff masses for the deflected grains is increased, so that more high-mass grains 
are prevented from entering the heliosphere and the m asses o f the dust p lumes would be 
enhanced compared to models based on compact grains. Hoang & Lazarianl (20141) find that 
larger grains are more sensitive to radiative torques, so that the polarizations of large huffy 
grains should be further enhanced over those of compact grains in the outer heliosheath. 


4. CONCLUDING REMARKS 

A filament of polarized starlight has been discovered within 10 pc of the heliosphere. The 
filament stars were initially identified from the rotation of polarization position angles with 
star distance. This study links these polarizations to the outer heliosheath regions where 
previous simulations indicate that interstellar nanograins are dehected in plume-like features 
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around the heliosphere. The filament polarizations show three primary observational prop¬ 
erties: (1) The best-fitting 1SMF vector to the filament polarization position angles coincides 
with the upwind velocity vector of interstellar neutrals flowing into the heliosphere. (2) The 
axis of this 100°-long filament is perpendicular to the BV plane, formed by the interstellar 
and magnetic and gas-velocity vectors, that shapes the heliosphere. (3) A contribution of 
±9.6° to the dispersion of the polarization position angles can be attributed to magnetic 
turbulence. 


We propose that the filament is formed in the outer heliosheath, where the 1SMF and 
interstellar plasma are deflected around the heliosphere. The polarizations would be found 
where the flow streamlines and draped 1SMF are approximately parallel to each other, and 
roughly perpendicular to the BV plane and sightline. However, although the filament traces 
several physical properties of the outer heliosheath, the difference in the thickness of the 
heliosheath compared to normal interstellar sightlines suggests that grain alignment mecha¬ 
nisms in the outer heliosheath would have to be more efficient than in the general ISM for 
such an origin. 

An outer heliosphere origin for the filament can be tested by evaluating the propaga¬ 
tion of porous interstellar dust grains through MHD models of the dynamically evolving 
outer heliosheath plasma and magnetic fields. Alignment efficiency in the outer heliosheath, 
compared to the ISM, would be augmented by increased grain charging, radiative torques, 
and stronger magnetic field strengths. Porous interstellar grains will substantially increase 
the mass and number of deflected interstellar grains above values for compact grains. The 
models will need to keep track of magnetic, collisional, and radiative torques on the grains in 
order to evaluate grain alignment and the asymmetric grain opacity that causes polarization. 
Filament properties that need to be explained include the filament width, the connection 
between the filament ISMF vector and neutral He flow vector, and detailed modeling of grain 
alignment mechanisms and polarization efficiencies for porous grains. 


The properties of the outer heliosheath are poorly known. While IBEX and the two 
Voyager missions have provided essential data on the outer heliosphere, data on the vector 
motion of the outer heliosheath plasma are not presently available because the Voyager 1 
plasma instrume nt i s not op erational and Voyager 2 has not yet reached the heliopause (e.g. 


Ri chardson & Decker! 12014). If an outer heliosheath origin for the filament polarizations is 


confirmed by theoretical modeling, a new diagnostic of the outer heliosheath will become 
available. This diagnostic would have the unique property of allowing the outer heliosheath 
to be monitored using ground-based observations. This convergence between interstellar 
dust data and heliosheath science would be useful for understanding the interactions be¬ 
tween astrospheres around external stars and the ISM, heliosphere-related foregrounds to 
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the cosmic microwave background, and the time-variable flux of interstellar grains into the 
inner heliosphere and terrestrial atmosphere. 


This research has been partly supported by the NASA Explorer program through sup¬ 
port for the IBEX mission, and by the European Research Council Advanced Grant HotMol 
(ERC-2011-AdG 291659). 
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Table 1. Stars tracing magnetic “filament” 


Star 

Coordinates 
b (deg) 

Distance 

(pc) 

Spec. 

Polarization A 
(10~ 5 ) 

$gal 

(deg) 

Ref. B 

Data used in 

Paper Ilf' 








HD 131977 


338, 

34 

6 

K4 V 

55.4 ± 24.0 

99.2 ±12.3 

LNA 

HD 161797 


52, 

26 

8 

G5 IV 

0.93 ±0.21 

92.3 ±6.7 

PP 

HD 172167 

a Lyr 

67, 

19 

8 

A0 Va 

1.72 ±0.1 

104.7 ± 1.4 

PP, W 

HD 159561 

a Oph 

36, 

23 

14 

A5 III V 

2.34 ±0.2 

96.4 ± 2.4 

PP, W 

HD 120467 


320, 

38 

14 

K1 III 

77.0 ± 23.0 

120.9 ±8.6 

San 

HIP 82283 


4, 

18 

18 

G6 III 

113.0 ±25.0 

100.1 ±6.3 

San 

HD 144253 


353, 

23 

19 

K3 V 

13.0 ±3.0 

115.2 ± 7.5 

NOT 

HD 119756 


316, 

29 

19 

F2 V 

90.0 ±35.0 

122.4 ± 10.6 

His 

HD 130819 


340, 

38 

24 

F3 V 

6.0 ±3.0 

114.9 ± 14.5 

NOT 

HD 134987 


339, 

27 

26 

G6 IV-V 

14.0 ±2.0 

100.6 ±4.5 

NOT 

HD 136894 


340, 

24 

28 

G8 V 

14.0 ±4.0 

106.6 ±8.0 

NOT 

HD 161096 


29, 

17 

25 

K2 III 

3.2 ±0.2 

90.9 ±1.9 

PP, W 

HD 161868 

7 Oph 

28, 

15 

29 

Al V 

4.1 ±0.3 

91.1 ±2.1 

PP, W 

New data from DiPol-2 

at UH88 

D 






HD 153631 


7, 

17 

26 

GO V 

9.0 ±0.9 

59.8 ±3.1 

UH e 

HD 141937 


352, 

27 

33 

G2 V 

7.9 ±1.0 

137.0 ±4.3 

UH 

HD 145518 


356, 

24 

33 

GO.5 V 

5.2 ± 1.1 

5.2 ± 1.1 

UH 


mhe units are degree of polarization. 


^References: LNA Data collected at LNA Picos dos Dios and listed in IFrisch et alJ (120120 or IFrisch c t al 


(2015 ); W-Wiktor owicz data from Lick observatory (2014, private co mmunication) ; PP PlanetPol data from 


Bailev et al. ( 2010l) ; NO T-Paper II; S an-LNA Picos dos Dios data from Santos et al. ( 2011 ); Hls-Compiled data 


from different sources in Heiled ( 2000l) .; UH-Data collected with the UH88 and presented in this paper 
c Frisch et al. ( 20151) 


D These data were obtained with DiPol-2 on the University of Hawaii 2.2m telescope in June and July 2014. 

e HD 153631 was also observed by NOT in 2010 ( Frisch et al.ll2012 ). giving P= 5±3 x 10 5 degree of polarization 
and 0 ga i= 61.2° ± 13.40° 







































22 


0.01 /zm grains 




Fig. 1.— The filament star polarizations (green bars) are shown plotted against the modeled 
column densities of radius 0.01 /jin interstellar dust grains interacting with the heliosphere. 
The small and large blue crosses indicate the direction of the IBEX ribbon ISMF and 
inflowing interstellar gas. The horizontal axis coincides with the ecliptic plane and the 
image is centered on the heliosphere nose defined by the inflowing ISM. The color coding 
shows the column density for 0.01 dust micron grains out to a distance of 400 AU. Color 
coding is based on the true column density (cm -2 ) for the grains for a gas-to-dust ratio of 
150. The dust distribution is from the dust trajectory models of ISlavin et al.l (1201 2) 
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Fig. 2.— The polarizations of the filament stars are compared to markers of the outer 
heliosheath region or beyond. Those polarizations, and the best-fitting magnetic field di¬ 
rection to those polarizations -Bfiiament, are plotted in galactic coordinates as black bars and 
the dotted-line box (which outlines the uncertainty range). The uncertainties on Filament 
overlap the upwind direction of the primary interstellar He° atoms (blue triangle) and the 
uncertainty range of the inflowing warm He° breeze (purple box, solid line). The ISMF 
di rections from the centers of the IBEX ribbon arc for the five energy bands of IBEX -Hi 
(IFunsten et al. 2013) are plotted as red open circles, with the size of the circle increasing with 
energy. The black diamonds r epresent t he locations of the low frequency outer heliosheath 
2-3 kHz emissions reported in iKurth fe Gurnettl (120031). The magnetic field direction that 
dominates polarization data (B CLI c, large black cross, Fri sch et al. 20151 ) overlaps the ISMF 
found from the IBEX ribbon. The dashed and solid gray lines show the locations of the 
BVplane and the axis of the filament, respectively. The lines enclose an angle of 80° ± 14° 
and are perpendicular to each other. The fluxes of the IBEX ribbon for the ~ 1.1 keV 
energy band are plotted on the color scale, with the hig hest fluxe s plotted as red and the 
lowest fluxes plotted in green (see Fig. 4 of ISchwadron et al.l 12014 1) . 
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Fig. 3.— Left: A cartoon showing the BV plane and filament star line plotted against a 
frontal view of the heliopause, i.e. as viewed from the upwind direction. The heli opaus e 


cartoon is borrowed from the MHD model displayed in Fig. la of flPo gorelov et al.l 2 0111). 


The BV plane corresponds to the hydrogen deflection plane displayed by flPogorelov et ah 


20 111) in Fig. la. The orientation of the filament star line is displayed by the draped lines. 


The filament orientation is perpendicular to the BV plane, and the filament magnetic field 
is directed toward the upwind velocity of interstellar neutrals, suggesting that a heliospheric 
origin requires the filament to form where the interstellar dust and magnetic field are fully 
deflected at the heliosphere nose. Right: The modeled trajectory of an interstellar grain 
such as may generate the polarization is displayed on the right, for a compact dust grain 
with radius 0.01 /irn approaching the heliosphere nose from the upwind direction. The 
grain gyroradius varies as the grain charge is modified by interactions with the heliospheric 
radiation and plasma components, until the grain is finally deflected ar ound the heliopause 
along with the interstellar plasma and magnetic field 
details). 


see 


Slavin et al.l 2012, for modeling 































